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ABSTRACT 

We study here schedulers for a class of rules that naturally 
arise in the context of rule-based constraint programming. 
We systematically derive a scheduler for them from a generic 
iteration algorithm of Apt Q. We apply this study to so- 
called membership rules of Apt and Monfroy . This leads 
to an implementation that yields for these rules a consid¬ 
erably better performance than their execution as standard 
CHR rules. 
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1. INTRODUCTION 

In this paper we are concerned with schedulers for a class of 
rules that naturally arise in the context of constraint pro¬ 
gramming represented by means of rule-based programming. 
An example of such rules are so-called membership rules, in¬ 
troduced in Apt and Monfroy j^. Their relevance stems 
from the following observations there made for constraint 
satisfaction problems (CSP’s) with finite domains: 

• constraint propagation can be naturally achieved by 
repeated application of the membership rules; 

• in particular the notion of hyper-arc consistency can 
be characterized in terms of the membership rules; 

• for constraints explicitly defined on small Hnite do¬ 
mains all valid membership rules can be automatically 
generated (For a most recent reference on the subject 
of such an automatic rule generation see Abdennadher 
and Rigotti Q.); 

• many rules of the CHR language (Constraint Handling 
Rules) of Friihwirth Q that are used in specific con¬ 
straint solvers are in fact membership rules. Now, 
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in the logic programming approach to constraint pro¬ 
gramming CHR is the language of choice to write con¬ 
straint solvers. 

In the resulting approach to constraint programming the 
computation process is limited to a repeated application of 
the rules intertwined with splitting (labeling). So the viabil¬ 
ity of this approach crucially depends on the availability of 
efficient schedulers for such rules. This motivates the work 
here reported. We provide an abstract framework for such 
schedulers and use it as a basis for an implementation. 

The abstract framework is based on an appropriate modi¬ 
fication of the generic approach to constraint propagation 
algorithms introduced in Apt Q and Apt Q . In this frame¬ 
work one proceeds in two steps. First, a generic iteration 
algorithm on partial orderings is introduced and proved cor¬ 
rect in an abstract setting. Then it is instantiated with spe- 
cffic partial orderings and functions to obtain specihc con¬ 
straint propagation algorithms. In this paper, as in Apt Q, 
we take into account information about the scheduled func¬ 
tions. Here we consider functions in the form of the rules 
b—>g, where b and g satisfy a number of natural conditions. 
We call such functions good rules. The relevant observation 
is that membership rules are good rules. Then we propose a 
specific scheduler in the form of an algorithm R, appropriate 
for good rules. 

The implementation is provided as an ECL*PS® program 
that accepts as input a set of membership rules and con¬ 
structs an ECL'PS® program that is the instantiation of the 
R algorithm for this set of rules. As membership rules can 
be naturally represented as CHR propagation rules, one can 
assess this implementation by comparing it with the perfor¬ 
mance of the standard implementation of membership rules 
in the CHR language. We found by means of various bench¬ 
marks that our implementation is considerably faster than 
CHR. 

CHR is available in a number of languages including the 
ECL'PS® and the Sicstus Prolog systems. In both cases 
CHR programs are compiled into the source language. There 
is also a recent implementation in Java, see Q. A great deal 
of effort was spent on implementing CHR efficiently. For an 
account of the most recent implementation see Holzbaur et 
al. 1^. Since, as already mentioned above, many CHR rules 
are membership rules, our approach provides a better im¬ 
plementation of a subset of CHR. While being stricly smaller 



than full CHR, the actual class of relevant rules is wider than 
the class of membership rules. The essential properties, such 
as monotonicity of condition and conclusion, are enjoyed 
by many rules that describe constraint propagation. This, 
hopefully, may lead to new insights into design and imple¬ 
mentation of languages appropriate for writing constraint 
solvers. 

It is important to stress that the discussed implementation 
was obtained by starting from “first principles” in the form 
of a generic iteration algorithm on an arbitrary partial or¬ 
dering. This shows the practical benehts of studying the 
constraint propagation process on an abstract level. 

2. REVISIONS OF THE GENERIC ITERA¬ 
TION ALGORITHM 
2.1 The Original Algorithm 

Let us begin our presentation with recalling the generic algo¬ 
rithm of Apt 1^. We slightly adjust the presentation to our 
purposes by assuming that the considered partial ordering 
also has the greatest element T. 

So we consider a partial ordering {D, C ) with the least 
element T and the greatest element T, and a set of functions 
F ~ {fi, ■ ■ ■, fk} on D. We are interested in functions that 
satisfy the following two properties. 

Definition 1. 

• / is called inflationary if x C f{x) for all x. 

• / is called monotonic x F y implies f(x) F f(y) 

for all X, y. □ 

Then the following algorithm is used to compute the least 
common hxpoint of the functions from F. 

d := T; 

G~F- 

while G yf 0 and d yf T do 

choose g £ G; 

G G — {g}; 

G G U update{G, g, d); 

d — g{d) 

end 

Figure 1: Generic Iteration Algorithm (GI) 

where for all G, g, d the set of functions update(G, g, d) from 
F is such that 

A {/ G F - G I /(d) = d A /(g(d)) yf g(d)} C update{G,g,d), 
B g{d) = d implies that update{G, g,d) = 0, 

C g{g{d)) yf g(d) implies that g £ update{G,g,d). 

Intuitively, assumption A states that update {G, g, d) con¬ 
tains at least all the functions from F — G for which the 
“old value”, d, is a hxpoint but the “new value”, g(d), is 
not. So at each loop iteration such functions are added to 


the set G. In turn, assumption B states that no functions 
are added to G in case the value of d did not change. As¬ 
sumption C provides information when g is to be added 
back to G as this information is not provided by A. On the 
whole, the idea is to keep in G at least all functions / for 
which the current value of d is not a hxpoint. 

The use of the condition d yf T, absent in the original pre¬ 
sentation, allows us to leave the while loop earlier. Our 
interest in the GI algorithm is clarihed by the following re¬ 
sult. 

Theorem 1 (Correctness). Suppose that all func¬ 
tions in F are inflationary and monotonie and that {D, C ) 
is finite and has the least element T and the greatest ele¬ 
ment T. Then every execution of the GI algorithm termi¬ 
nates and computes in d the least common fixpoint of the 
functions from F. 

Proof. (Sketch). The following statement is an invariant 
of the while loop of the algorithm: 

(yf^F-G f{d) = d) A (V/ e F /(T) = T). 

This implies that the algorithm computes in d a common 
hxpoint of the functions from F. The fact that this is the 
least common hxpoint follows from the assumption that all 
functions are monotonic. 

In turn, termination is established by considering the lexi¬ 
cographic ordering of the strict partial orderings (!),□) and 
(Af, <), dehned on the elements of F x A/” by 

(di,ni) <ie 3 , (d 2 ,n 2 ) ih di □ d 2 or (di = d 2 and ni < n 2 ). 

Then with each while loop iteration of the algorithm the 
pair {d,card G), where card G is the cardinality of the set 
G, strictly decreases in the ordering <iex- D 

2.2 Removing Functions 

We now revise the GI algorithm by modifying dynamically 
the set of functions that are being scheduled. The idea is 
that, whenever possible, we remove functions from the set 
F. This will allow us to exit the loop earlier which speeds 
up the execution of the algorithm. 

To realize this idea we proceed as follows. First, we in¬ 
troduce the following property that will be satished by the 
considered functions. 

Definition 2. Suppose d £ D and f £ F. We say that / 
is stable above d if d □ e implies /(e) = e. We then say that 
/ is stable if it is stable above /(d), for all d. □ 

That is, / is stable if for all d and e, /(d) C e implies 
/(e) = e. So stability implies idempotence, which means 
that /(/(d)) = /(d), for all d. Moreover, if d and /(d) are 
comparable for all d, then stability implies inflationarity. In¬ 
deed, if d C /(d), then the claim holds vacuously. And if 
/(d) C d, then by stability /(d) = d. 

Next, we assume that for each function g £ F and each 
element d £ D, two lists of functions from F are given. 


friends{g,d) and obviated{g,d) that satisfy the following 
condition 

VdVe □ g ogi o ... o g^{d) 

V/ € friends{g,d) U obviated{g,d) (/(e) = e) (1) 
where friends{g, d) = [gi,..gt]- 

That is, for all d, each function / in friends{g,d) U 
obviated(g, d) is stable above g o gi o ... o gk(d). 

Now, we modify the GI algorithm in such a way that each ap¬ 
plication of 5 to d will be immediately followed by the appli¬ 
cations of all functions from friends{g,d) and by a removal 
of the functions from friends{g, d) and from obviated{g,d) 
both from F and G. This modified algorithm is shown in 
Fig. To keep the notation uniform we identified at some 
places the lists friends{g, d) and obviated{g, d) with the sets. 


2.3 Functions in the Form of Rules 

In what follows we consider the situation when the scheduled 
functions are of a specific form b g, where 6 is a condition 
and g a function, that we call a body. We call such functions 
rules. 

First, we explain how rules are applied. Given an element 
d of D, a condition b evaluates in d to either true or false, 
denoted Holds{b, d) and -^Holds{b,d), resp. 

Given a rule b g we define then its application as follows: 

— / 3{d) if Holds{b,d) 

— ii^Holdsib,d) . 

The rules introduced in the next section will be of a specific 
type. 


d := T; 

Fo ~ F-, 

G~F- 

while G 7 ^ 0 and d 7 ^ T do 

choose g £ G; 

G G — {g}; 

F ~ F — {friends{g,d) U obviated{g,d)); 

G ~ G — {friends{g, d) U obviated{g, d)); 

G ~ G U update{G, h, d), 

where h = g o gi o ... o g^ and friends{g, d) = [gi ,..., gt]; 
d := h(d) 

end 

Figure 2: Revised Generic Iteration Algorithm (RGI) 

The following result then shows correctness of this algo¬ 
rithm. 

Theorem 2. Suppose that all functions in F are infla¬ 
tionary and monotonic and that {D, C ) is finite and has the 
least element T and the greatest element T. Additionally, 
suppose that for each function g £ F and d £ D two lists of 
functions from F are given, friends{g,d) and obviated{g, d) 
such that condition (nW holds. 

Then the Correctness Theorem holds with the GI algorithm 
replaced by the RGI algorithm. 

Proof. In view of condition (|^) the following statement 
is an invariant of the while loop: 


Definition 3. Consider a partial ordering {D, C ). 

• We say that a condition b is monotonic if Holds{b, d) 
and d G e implies Holds{b,e), for all d, e. 

• We say that a condition b is precise if the least d exists 
such that Holds(b,d). We call then d the witness for 
b. 

• We call a rule b ^ g good if h is monotonic and precise 

and g is stable. □ 

When all rules are good, we can modify the RGI algorithm by 
taking into account that an application of a rule is a two step 
process: testing of the condition followed by a conditional 
application of the body. This will allow us to construct the 
lists friends{g, d) and obviated{g,d) before the execution of 
the algorithm, without using the parameter d. Moreover, 
the list friends{g) can be constructed in such a way that 
the conditions of its rules do not need to evaluated at the 
moment they are applied, as they will all hold. The details 
of a specific construction that we shall use here will be given 
in a moment, once we identify the condition that is crucial 
for the correctness. This revision of the RGI algorithm is 
given in Fig. |^. 

Again, we are interested in identifying conditions under 
which the Correctness Theorem holds with the GI algo¬ 
rithm replaced by the R algorithm. To this end, given a rule 
fe —» 3 in F and d £ D, define as follows: 


yf£F-G (fid) = d) A V/ e F (/(T) = T) A 

V/ e Fo - F Ve □ d (/(e) = e). (2) 


friends{b —> g, d) 


friends{b —> g) if Holds{b, d) 
[] a ^Holds{b,d) 


So upon termination of the algorithm the conjunction of this 
invariant with the negation of the loop condition, i.e., 

G = 0 V d = T 

holds, which implies that V/ £ Fo (/(d) = d). 

The rest of the proof is the same. □ 


obviated {b 


g,d) := 


obviated{b 

[b^g] 


g) if Holds{b, d) 
if Ve □ d 
-^Holds{b, e) 
otherwise 


We now have the following counterpart of the Gorrectness 
Theorem 


In the next section we shall focus on functions that are in 
a special form. For these functions we shall show how to 
construct specific lists friends{g, d) and obviated(g, d). 


Theorem 3 (Correctness). Suppose that all func¬ 
tions in F are good rules of the form b ^ g, where g is infla¬ 
tionary and monotonic, and that {D, C ) is finite and has 



d := _L; 

Fo := F; 

G~F- 

while G 7^ 0 and d 7^ T do 

choose f £ G', suppose f is b ^ g; 

G~G-{b^g}-, 
if Holds {b, d) then 

F ~ F — {friends{b g)U obviated(b ^ g))\ 

G ~ G — {friends{b —> g) U obviated{b g)); 

G ~ G U update{G, h, d), 
where h = g o g^ o . . . o g^ 
and friends{b ^ g) = [bi ^ gi,.. .,bk ^ fffc]; 
d ■- h{d) 
else 

if Ve □ d -iHolds{b,e) then 
F~F-{b^g} 

end 

end 

end 

Figure 3: Rules Algorithm (R) 

the least element _L and the greatest element T. Further, as¬ 
sume that for each rule b—>g the lists friendsib —> g, d) and 
obviated(b ^ g,d) defined as above satisfy condition and 
the following condition: 

\/d{b'^ g'£ friends {b ^ g) A ( 3 ) 

Holds{b, d)^Ve^ g{d) Holdsib', e) . 

Then the Correctness Theorem holds with the GI algorithm 
replaced by the R algorithm. 

Proof. It suffices to show that the R algorithm is an in¬ 
stance of the RGI algorithm. On the account of condition (^) 
and the fact that the rule bodies are inflationary functions, 
Holds{b, d) implies that 

((b^g)o(bi ^ gi)o. . .o(bk gk))id) = (gogio. . .ogk){d), 

where friends{b g) = [61 —> yi,..., ^ gfc]. This takes 

care of the situation when if Holdsib, d). 

In turn, the definition of friendsfi ^ g,d) and 
obviated(b ^ g,d) and assumption B take care of the 
situation when if -^Holds{b,d). When the condition b fails 
for all e □ d, then we can conclude that for all such e we 
have {b—fg)ie) = e. This allows us to remove at that 
point of the execution the rule b ^ g from the set F. This 
amounts to adding 6 —> 5 to the set obviated(b g,d) at 
runtime. Note that condition (Q) is then satisfied. □ 

We now provide an explicit construction of the lists friends 
and obviated for a rule b g in the form of the algorithm 
in Fig. ^ Gl(d) stands here for the GI algorithm activated 
with T replaced by d and the considered set of rules as the 
set of functions F. Further, given an execution of Gl(e), we 
call here a rule g relevant if at some point g(d) 7^ d holds 
after the “choose g £ G” action. 

Note that b^g 0 friendsib ^ g) since b^g is a good 
rule, while b g £ obviatedib —> g) since by the stability of g 
(/(e) = e holds. 


e ;= witness of 6; 
e := Gl(g(e)); 

friendsib ^ g) ~ list of the relevant rules h £ F 
in the execution of Gl((/(e)); 
obviatedib —>(?):=[]; 

for each (F g') £ F — friends —> g) do 
if g'ie) = e or Ve' □ e -^Holdsib', e') then 
obviatedib g) := \b' —> g'\obviatedib (/)] 

end 

end 

Figure 4: Friends and Obviated Algorithm (F & 0) 

The following observation now shows the adequacy of the 
F & D algorithm for our purposes. 

Lemma 1. Uwn termination of the F & 0 algorithm con¬ 
ditions and hold, where the lists friendsib g, d) and 
obviatedib ^ g,d) are defined as before Theorem^. □ 

Let us summarize now the findings of this section that cul¬ 
minated in the R algorithm. Assume that all functions are 
of the form of the rules satisfying the conditions of the Cor¬ 
rectness Theorem Then in the R algorithm, each time 
the evaluation of the condition b of the selected rule b g 
succeeds, 

• the rules in the list friendsib g) are applied directly 
without testing the value of their conditions, 

• the rules in friendsib —> g) U obviatedib —> g) are per¬ 
manently removed from the current set of functions G 
and from F. 

2.4 Recomputing of the Least Fixpoints 

Another important optimization takes place when the R al¬ 
gorithm is repeatedly applied to compute the least fixpoint. 
More specifically, consider the following sequence of actions: 

• we compute the least common fixpoint d of the func¬ 
tions from F, 

• we move from d to an element e such that d C e, 

• we compute the least common fixpoint above e of the 
functions from F. 

Such a sequence of actions typically arises in the framework 
of CSP’s, further studied in Section ^ The computation of 
the least common fixpoint d of the functions from F cor¬ 
responds there to the constraint propagation process for a 
constraint G. The moving from d to e such that d C e corre¬ 
sponds to splitting or constraint propagation involving an¬ 
other constraint, and the computation of the least common 
fixpoint above e of the functions from F corresponds to an¬ 
other round of constraint propagation for G. 

Suppose now that we computed the least common fixpoint 
d of the functions from F using the RGI algorithm or its 
modification R for the rules. During its execution we per¬ 
manently removed some functions from the set F. Then 



these functions are not needed for computing the least com¬ 
mon hxpoint above e of the functions from F. The precise 
statement is provided in the following simple, yet crucial, 
theorem. 

Theorem 4. Suppose that all functions in F are infla¬ 
tionary and monotonic and that {D, F ) is finite. Suppose 
that the least common fixpoint do of the functions from F is 
computed by means of the RGI algorithm or the R algorithm. 
Let Ffin be the final value of the variable F upon termination 
of the RGI algorithm or of the R algorithm. 

Suppose now that do F e. Then the least common fixpoint 
eo above e of the functions from F coincides with the least 
common fixpoint above e of the functions from Fg^. 

Proof. Take a common fixpoint ei of the functions from 
Ffi„ such that e O ei. It suffices to prove that ei is common 
fixpoint of the functions from F. So take f £ F — Ffin. 
Since condition (0) is an invariant of the main while loop 
of the RGI algoritnm and of the R algorithm, it holds upon 
termination and consequently / is stable above do. But 
do F e and e C ei, so we conclude that /(ei) = ei. □ 

Intuitively, this result means that if after splitting we re¬ 
launch the same constraint propagation process we can dis¬ 
regard the removed functions. 

In the next section we instantiate the R algorithm by a set of 
rules that naturally arise in the context of constraint satis¬ 
faction problems with Hnite domains. In Section^ we assess 
the practical impact of the discussed optimizations. 

3. CONCRETE FRAMEWORK 

We now proceed with the main topic of this paper, the sched¬ 
ulers for the rules that naturally arise in the context of con¬ 
straint satisfaction problems. First we recall briefly the nec¬ 
essary background on the constraint satisfaction problems. 

3.1 Constraint Satisfaction Problems 

Consider a sequence of variables X := xi,...,Xn where 
n > 0, with respective domains Di,..., Dn associated with 
them. So each variable Xi ranges over the domain Di. By 
a constraint C on X we mean a subset of Di x ... x Dn. 
Given an element d := di ,.. .,dn of Di x ... x and a 
subsequence Y := of X we denote by d[Y] the 

sequence di,^, .. .,dig. In particular, for a variable Xi from 
X, d\xi\ denotes di. 

Recall that a constraint satisfaction problem, in short CSP, 
consists of a finite sequence of variables X with respec¬ 
tive domains D, together with a Hnite set C of con¬ 
straints, each on a subsequence of X. We write it as 
{C ■, xi € Di,.. .,x„ £ Dn), where X ~ xi,...,Xn and 
V := Di,..., Dn. 

By a solution to {C ■, xi £ Di,...,Xn £ Dn) we mean an 
element d £ D\ x ... x Dn such that for each constraint 
C £ C on a sequence of variables X we have d[X] £ C. We 
call a CSP (insistent if it has a solution. Two CSP’s with 
the same sequence of variables are called equivalent if they 
have the same set of solutions. 


3.2 Partial Orderings 

With each CSP V '.= {C xi £ Di,...,Xn £ Dn) 
we associate now a specific partial ordering. Initially 
we take the Cartesian product of the partial orderings 
{V{Di), D),..., {V{Dn), 3). So this ordering is of the form 

(P(T>i) X ... X V{Dn),^) 

where we interpret D as the the Cartesian product of the 
reversed subset ordering. The elements of this partial or¬ 
dering are sequences {Ei,..., En) of respective subsets of 
(Di,..., Dn) ordered by the componentwise reversed subset 
ordering. Note that in this ordering (Hi,..., Dn) is the least 
element while 

n times 

is the greatest element. However, we would like to iden¬ 
tify with the greatest element all sequences that contain as 
an element the empty set. So we divide the above partial 
ordering by the equivalence relation i ?0 according to which 

{El,...,En) Rfl){Fl,...,Fn) iff 
(Hi,.. .,En) = (Hi,..., En) or (3i Ei = 9 and 3j Fj = 0). 

It is straightforward to see that i ?0 is indeed an equivalence 
relation. 

In the resulting quotient ordering there are two types of 
elements: the sequences (Hi,...,H„) that do not contain 
the empty set as an element, that we continue to present in 
the usual way with the understanding that now each of the 
listed sets is non-empty, and one “special” element equal to 
the equivalence class consisting of all sequences that contain 
the empty set as an element. This equivalence class is the 
greatest element in the resulting ordering, so we denote it 
by T. In what follows we denote this partial ordering by 
[Dr, C ). 

3.3 Membership Rules 

Fix now a specific CSP V {C ; xi £ Hi ,.. .,Xn £ Dn) 
with finite domains. We now recall the rules introduced in 
Apt and Monfroy They are called membership rules 

and are of the form 

y\ £ Si,.. .,yk £ Sk ^ zi ^ ai , . . Zm 

where 

• yi,.. .,yk are pairwise different variables from the set 
{xi,..., Xn} and Si,..., Sj, are subsets of the respec¬ 
tive variable domains, 

• zi,.. .,Zm are variables from the set { 2 : 1 ,..., Xn} and 
tti,.. .,am are elements of the respective variable do¬ 
mains. 

Note that we do not assume that the variables zi,.. .,Zm are 
pairwise different. 

The computational interpretation of such a rule is: 

^In our presentation we slightly relax the original syntactic 
restrictions. 



if for i £ the current domain of the variable 
Ui is included in the set Si, then for j £ 
remove the element ai from the domain of Zi. 


is a membership rule and when for each rule b ^ g the lists 
friends{b—f g) and obviated{b ^ g) are constructed by the 
F & 0 algorithm. 


When each set Si is a singleton, we call a membership rule 
an equality rule. 


Note 1 . Consider the partial ordering {D-p, C ). 


Let us reformulate this interpretation so that it hts the 
framework considered in the previous section. To this end 
we need to clarify how to evaluate a condition, and how to 
interpret a conclusion. We start with the first item. 


Definition 4 - Given a variable y with the domain Dy and 
E,S D Dy we dehne 


Holdsiy €S,E) iff EC S, 


and extend the definition to the elements of the considered 
ordering {Dp, C ) by putting 


(i) Each membership rule is good. 

(ii) Each function zi 7^ ai, .. .,Zm 7^ Um on Dp is 

• inflationary, 

• monotonic. □ 


To be able to instantiate the algorithm R with the member¬ 
ship rules we still need to define the set update{G,g,d). In 
our implementation we chose the following simple dehnition: 


update{G, b—fg,d) 


F -G if Holds{b, d) A g{d) / d 
0 otherwise. 


Holds{y €S,{Ei,.. ., Er,)) iff Ek C S', 
where we assumed that y is Xk, 

Holdsiy e *S',T). 

Then we interpret a sequence j/i £ Si,..., yj, £ Sj, of condi¬ 
tions as a conjunction, so by putting 

Holds{yi € Si,.. .,yk £ Sk, {Ei,.. .,En)) iff 
Holdsiyi £ Si, {El,.. .,En)) for i £ [l..fc], 

and 

Holds{yi € Si,...,yk C Sk,T). □ 

Concerning the second item we proceed as follows. 

Definition 5 . Given a variable 2: with the domain Dz we 
interpret the atomic formula 2 7^ a as a function on ViDz), 
dehned by: 

{z^a){E) :=£-{«}. 

Then we extend this function to the elements of the consid¬ 
ered ordering {Dp, C ) as follows: 

• on the elements of the form {Ei,...,En) we use 
“padding”, that is we interpret it as the identity on 
the other components. If the resulting sequence con¬ 
tains the empty set, we replace it by T, 

• on the element T we put (2 7^ u)(T) := T 

Finally, we interpret a sequence 21 7^ ai,...,2„i 7^ Um of 
atomic formulas by interpreting each of them in turn. □ 

In view of the Gorrectness Theorem ^ the following observa¬ 
tion allows us to apply the R algorithm when each function 


To illustrate the intuition behind the use of the lists 
friends{b g) and obviated{b —> g) take the CSP V ~ 

(C ; xi £ {a,6, c},a:2 £ {a,b,c},X3 £ {a,b,c},X 4 £ {a, 6, c}) 

and consider the membership rules 

ri := xi £ {a, b} X2 a, xi 7^ b, 

r2 ■= xi £ {a, b}, X2 £ {b, c} ^ X3 a, 

rs := X 2 £ {b} —> *3 7 ^ 0,2:4 7^ b. 

Then upon application of rule ri rule r2 can be applied 
without evaluating its condition and subsequently rule rs 
can be deleted without applying it. So we can put rule r2 
into friends{ri) and rule rs into obviated{ri), and this in 
fact is what the F & 0 algorithm does. 

4. IMPLEMENTATION 

In this section we discuss the implementation of the R al¬ 
gorithm for the membership rules and compare it by means 
of various benchmarks with the CHR implementation in the 
ECL^PS'' system. 

4.1 Modelling of the Membership Rules in CHR 

Following Apt and Monfroy 0 the membership rules are 
represented as CHR propagation rules with one head. Recall 
that the latter ones are of the form 

Gi,...,Gi\ Bi,...,Bm. 

where 

• I > 0 , m > 0 , 

• the atom H of the head refers to the defined con¬ 
straints, 

• the atoms of the guard Gi,... ,Gi refer to Prolog re¬ 
lations or built-in constraints, 

• the atoms of the body Bi,..., Bm are arbitrary atoms. 



Further, recall that the CHR propagation rules with one head 
are executed as follows. First, given a query (that represents 
a CSP) the variables of the rule are renamed to avoid vari¬ 
able clashes. Then an attempt is made to match the head 
of the rule against the first atom of the query. If it is suc¬ 
cessful and the guard of the instantiated version of the rule 
succeeds, the instantiated version of the body of the rule is 
executed. Otherwise the next rule is tried. 

Finally, let us recall the representation of a membership rule 
as a CHR propagation rule used in Apt and Monfroy |^. Con¬ 
sider the membership rule 

yi € Si,.. .,yk e Sk ^ zi ^ ai , . . Zm 

related to the constraint c on the variables Xi,...,X„. 
We represent it as a CHR rule with the single head 
atom c{Xi,.. .,Xn) and guard atoms in{yi,Si) where 
the in/2 predicate is defined by in(X,L) dom(X,D), 
subset(D,L). The body consists of atomic calls Zi ##ai. 

In general^the application of a membership rule as defined 
in Section H and the execution of its representation as a CHR 
propagation rules coincide. Moreover, by the semantics of 
CHR, the CHR rules are repeatedly applied until a fixpoint is 
reached. So a repeated application of a finite set of member¬ 
ship rules coincides with the execution of the CHR program 
formed by the representations of these membership rules as 
propagation rules. 

4.2 Benchmarks 

In our approach the repeated application of a finite set of 
membership rules is realized by means of the R algorithm 
of Section g implemented in ECL*PS'^. The compiler con¬ 
sists of about 1500 lines of code. It accepts as input a set 
of membership rules, each represented as a CHR propagation 
rule, and constructs an ECL'PS'^ program that is the instan¬ 
tiation of the R algorithm for this set of rules. As in CHR, for 
each constraint the set of rules that refer to it is scheduled 
separately. 

For each considered constraint we use rules generated by a 
program discussed in 0. Our compiler constructs then for 
each rule g the lists friends{g) and obviated{g) by executing 
the F & 0 algorithm (essentially computing a fixpoint for 
each rule). Time spent on this construction is comparable 
with rule generation time. 

We chose benchmarks that embody several successive prop¬ 
agation steps, i.e., propagation interleaved with domain 
splitting or labelling. In Table ^ we list the results for se¬ 
lected single constraints. For each such constraint, say C on 
a sequence of variables xi,.. .,Xn with respective domains 
Di,.. Dn, we consider the CSP {C ■, xi £ Di, ...,*„ € Dn) 
together with randomized labelling. That is, the choices of 
a variable, value, and an assignment or a removal of the 
value, are random. The computation of only the solutions 
yields times that are insignificant, so the test program com¬ 
putes also all intermediate fixpoints, where some domains 
are not singleton sets. Branching at these recorded points 
takes place only once, that is, backtracking occurs immedi¬ 
ately if a recorded point is encountered again. In Table 
we report the results for CSP’s that formalize sequential au¬ 
tomatic test pattern generation for digital circuits (ATPG). 


These are rather large CSP’s that employ the and constraints 
of Table ^ and a number of other constraints. They are taken 
from a recent study by the first author that will be reported 
elsewhere. 

We measured the execution times for three rule schedulers: 
the standard CHR representation of the rules, the generic 
chaotic iteration algorithm GI, and its improved derivative R. 
The codes of the latter two algorithms are both produced 
by our compiler and are structurally equal, hence allow a 
direct assessment of the improvements embodied in R. 

An important point in the implementations is the question 
of when to remove solved constraints from the constraint 
store. The standard CHR representation of membership rules 
does so by containing, beside the propagation rules, one CHR 
simplification rule for each tuple in the constraint definition. 
Once its variables are assigned values that correspond to a 
tuple, the constraint is solved, and removed from the store 
by the corresponding simplification rule. This ‘solved’ test 
takes place interleaved with propagation. The implementa¬ 
tions of GI and R check after closure under the propagation 
rules. The constraint is considered solved if all its variables 
are fixed, or, in the case of R, if the set F of remaining rules 
is empty. 

In the tables we provide for each constraint or CSP the ratio 
of the execution times in seconds between, first, R and GI, 
and second, R and CHR. This is followed by the absolute times 
for R and GI / CHR. 


Const. 

rcc8 

fork 

and3 

and9 

andll 

MEM 

rel. 

26%/ll% 

43%/40% 

58%/47% 

13%/6% 

13%/3% 

abs. 

109 

0.23 

0.22 

70 

55.6 


419/950 

0.54/0.58 

0.38/0.47 

534/1096 

427/2077 

EQU 

rel. 

95%/100% 

95%/89% 

82%/74% 

94%/97% 

89%/94% 

abs. 

323 

18.9 

0.31 

286 

299 


341/324 

19.9/21.2 

0.38/0.42 

303/294 

335/318 


Table 1 : Randomized search trees for constraints 


Logic 

3-valued 

9-valued 

11-valued 7 

MEM 

relative 

64% / 35% 

71% / 24% 

85% / 86% 

absolute 

1.39 2.16/4.01 

124 175/509 

797 933/3120 

EQU 

relative 

63% / 70% 

44% / 59% 

39% / 48% 

absolute 

0.72 1.15/2.58 

2.40 5.50/4.09 

12.3 31.6/25.7 


Table 2 : CSP’s formalizing sequential ATPG 


4.3 Recomputing of the Least Fixpoints 

Finally, let us illustrate the impact of the permanent removal 
of the rules during the least fixpoint computation, achieved 
here by the use of the lists friends{g) and obviated{g). Given 
a set F of rules call a rule g £ F solving if friends (g) U 
obviated{g) = F. 

Take now as an example the equivalence relation = from 
three valued logic of Kleene j^[page 334] that consists of 












three values, t (true), f (false) and u (unknown). It is defined 
by the truth table 


= 

t 

f 

u 

t 

t 

f 

u 

f 

f 

t 

u 

u 

u 

u 

u 


The program of Apt and Monfroy generates for it 26 
minimal valid membership rules. Out of them 12 are solving 
rules. For the remaining rules the sizes of the set friends U 
obviated are; 17 (for 8 rules), 14 (for 4 rules), and 6 (for 2 
rules). 

In the R algorithm a selection of a solving rule leads directly 
to the termination {G = 0) and to a reduction of the set 
F to 0. For other rules also a considerable simplification 
in the computation takes place. For example, one of the 8 
rules with 17 rules in its set friends U obviated is 

r := X £ {0}, 2 € {0, m} —> 3 / / 0. 

Consider the CSP {= ; x £ {0}, y £ {0,1, w}, 2 £ {0, w}). In 
the R algorithm the selection of r is followed by the applica¬ 
tion of the rules in friends and the removal of the rules in 
friends U obviated. This brings the number of the considered 
rules down to 26 — 17 = 9. The R algorithm subsequently 
discovers that none of these nine rules is applicable at this 
point, so this set F remains upon termination. Then in 
a subsequent constraint propagation phase, launched after 
splitting or after constraint propagation involving another 
constraint, the hxpoint computation by means of the R al¬ 
gorithm involves only these nine rules instead of the initial 
set of 26 rules. For solving rules, this fixpoint computation 
immediately terminates. 

Interestingly, as Table ^ shows, the solving rules occur quite 
frequently. We list there for each constraint and each type of 
rules the number of solving rules divided by the total number 
of rules, followed in a new line by the average number of rules 
in the set friends{g) U obviated{g). 


Eind2 

andS 

and9 

andll 

fork 

rcc 8 

alien 

EQU 6/6 

13/16 113/134 

129/153 

9/12 183/183 498/498 

6 

14 

130 

148 

11 

183 

498 

MEM 6/6 

4/13 

72/1294 196/4656 0/24 

0/912 

-/26446 

6 

7 

810 

3156 

9 

556 

- 


Table 3 : Solving rules 


The fork constraint is taken from the Waltz language for the 
analysis of polyhedral scenes. The rcc 8 is the composition 
table for the Rerion Connection Calculus with 8 relations 
from Egenhofer |H. It is remarkable that all its 183 mini¬ 
mal valid equality rules are solving. While none of its 912 
minimal valid membership rule for rcc 8 is solving, on the 
average the set friends (g) U obviated (g) contains 556 mem¬ 
bership rules. Also all 498 minimal valid equality rules for 
the alien constraint, that represents the composition table 
for Allen’s qualitative temporal reasoning, are solving. 


The number of minimal valid membership rules exceeds 
26,000 and consequently they are too costly to analyze. 

The savings obtained by means of the lists friends (g) and 
obviated (g) are orthogonal to the ones obtained by a trans¬ 
formation of the CHR propagation rules into the sirnnlihca- 
tion rules discussed in Abdennadher and Rigotti [0]. We 
think that there is a relation between two approaches that 
we plan to study closer. 
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